Dynamical analogy between epileptic seizures and seismogenic 
electromagnetic emissions by means of nonextensive statistical mechanics 

Konstantinos Eftaxias*^, George Minadakis^ , Stelios. M. Potirakis'^, George Balasis'^ 

"Department of Physics, Section of Solid State Physics, University of Athens, Panepistimiopolis, GR 15784, Zografos, 

Athens, Greece 

''Department of Electronic and Gomputer Engineering, Brunei University, Kingston Lane, Uxbridge, Middlesex, UBS 3PH, 

U.K. 

"^Department of Electronics, Technological Educational Institute of Piraeus, 250 Thivon & P. Ralli, GR- 12244> Aigaleo, 

Athens, Greece 

'^Institute for Space Applications and Remote Sensing, National Observatory of Athens, Metaxa and Vas. Pavlou St., 

Penteli, 15236 Athens, Greece 



Abstract 

The field of study of complex systems considers that the dynamics of complex systems are founded on 
universal principles that may be used to describe a great variety of scientific and technological approaches 
of different types of natural, artificial, and social systems. Several authors have suggested that earthquake 
dynamics and neurodynamics can be analyzed within similar mathematical frameworks. Recently, authors 
have shown that a dynamical analogy supported by scale-free statistics exists between seizures and earth- 
quakes, analysing populations of different seizures and earthquakes, respectively. The purpose of this paper 
is to suggest a shift in emphasis from the large to the small scale: our analyses focus on a single epileptic 
seizure generation and the activation of a single fault (earthquake) and not on the statistics of sequences 
of different seizures and earthquakes. We apply the concepts of the nonextensive statistical physics to sup- 
port the suggestion that a dynamical analogy exists between the two different extreme events, seizures and 
earthquakes. Wc also investigate the existence of such an analogy by means of scale-free statistics (the 
Gutenberg-Richter distribution of event sizes and the distribution of the waiting time until the next event). 
The performed analysis confirms the existence of a dynamic analogy between earthquakes and seizures, 
which moreover follow the dynamics of magnetic storms and solar fiares. 

Keywords: nonextensivity, complex system dynamics, preseismic electromagnetic emissions, epileptic 
seizures 



1. Introduction 

Complexity nowadays is a frequently used yet poorly defined - at least quantitatively speaking - concept. 
It tries to embrace a great variety of scientific and technological approaches of all types of natural, artificial 
and social systems [1]. In this direction several examples have been given providing the elements of such 
definition [2-9] . When one considers a phenomenon that is "complex" refers to a system whose phenomeno- 
logical laws, which describe the global behaviour of the system, are not necessarily directly related to the 
"microscopic" laws that regulate the evolution of its elementary parts [10]. The field of study of complex 
systems considers that the dynamics of complex systems are founded on universal principles that may used 
to describe disparate problems [2]. This is a basic reason for our interest in complexity [6, 10-13]. There is 
an apparent paradox in the above mentioned suggestion. How is it possible for a concept as multifaceted 
as complexity to serve as a unifying direction? There is a common factor in these seemingly diverse phe- 
nomena. The complex systems adopt a pattern of behaviour almost completely determined by the collective 
effects. They exhibit remarkable properties of self-organization and emergence of coherent structure over 
many scales. The main feature of collective behaviour is that an individual unit's action is dominated by 
the influence of its neighbours; the unit behaves differently from the way it would behave on its own, so 
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that, all units simultaneously alter their behaviour to a common pattern. Thus, new features emerge as we 
move from one scale to another and the science of complexity is about revealing the principles that govern 
the ways in which new properties appear [10]. 

Complexity measures are generally performed to study pattern, structure and correlation of systems 
and processes. Although the definition of the term 'complexity' is not unique depending on the different 
contexts and different fields of application, its quantification has received considerable attention [14, and 
references therein]. In any case quantities such as: the disequilibrium, the randomness or uncertainty, the 
degree of order or organization, have played significant role on giving the notion of complexity [14, and 
references therein]. Complex systems seem to occur close, in some sense, to the frontier between order and 
disorder where most of their basic quantities exhibit nonexponential behaviours, very frequently power laws. 
It happens that the distributions and other relevant quantities that emerge naturally within the frame of 
noncxtensive statistical mechanics are precisely of this type [1]. Herein our analysis is mainly focused on 
the noncxtensive Tsallis statistical mechanics, namely, a recently introduced noncxtensive model for EQ 
dynamics [15] and Tsallis entropy. One universal footprint seen in many complex systems is self-affinity and 
the fractional power law relationship is a classic expression of a self-affine structure - the fractal. Note that 
the noncxtensive model also leads to a fractional power law relationship [16, 17]. 

The contribution in this work enhances the suggestion that transferring ideas, methods and insights from 
investigations in hitherto disparate areas will cross-fertilize and lead to important new results. We stress 
that time-clustering methods are so powerful in describing so diverse natural as well as human phenomena. 
For instance, natural phenomena as well as human activities were found to share common time-clustering 
behaviour [18-24]. 

Epileptic seizures (ESs) and earthquakes (EQs) are complex phenomena, which have highly intricate 
cluster and hierarchical structures, spatial and temporal correlation with feedback, self-organization and 
connection diversity. Authors have suggested that dynamics of EQs and neurodynamics can be analyzed 
within similar mathematical frameworks [25, 26]. Characteristically, driven systems of interconnected blocks 
with stick-slip friction capture the main features of EQ process. These models, in addition to simulating 
the aspects of EQs and frictional sliding, may also represent the dynamics of neurological networks [25]. 
Hopfield [27] proposed a model for a network of N integrate-and-fire neurons. In this model, the dynamical 
equation of /c*'* neuron, see equation 28 in [27] is based on the Hodgekin-Huxley model for neurodynamics 
and represents the same kind of mean field limit that has been examined in connection with EQs [26]. 

Recently, Osorio et al. [28] in a pioneering work have shown that a dynamical analogy supported by 
scale-free statistics exists between ESs and EQs. More precisely, the authors performed the analysis using: 
(i) a population of different EQs between 1984-2000 available in the Southern California Seismic Network 
catalogue, and (ii) a population of different ESs from continuous multiday voltage recordings directly from 
the brains of 60 human subjects with epilepsy undergoing surgical evaluation at the University of Kansas 
Medical Center between 1996 and 2000. In this work we investigate the existence of a dynamical analogy 
between ESs and EQs at the level of a single fault / seizure activation, namely, we examine whether a 
dynamical analogy may exist for the ways in which firing neurons / opening cracks organize themselves to 
produce a single ES / EQ. 

We perform the analysis mainly based on a recently introduced noncxtensive model for EQ dynam- 
ics which leads to a Gutenberg-Richter type law for the relationship between frequency and magnitude 
of EQs [15, 29]. We show that the populations of: (i) electric pulses included in a single ES, (ii) fracto- 
electromagnctic pulses rooted in the activation of a single fault, (iii) different ESs occurred in different 
human's brains, (iv) different EQs occurred in different faults in various seismic regions, follow the above 
mentioned noncxtensive statistical law. The similarity is extended up to the laboratory seismicity. Impor- 
tantly, the noncxtensive model is also able to describe solar flares and magnetic storms, as well. 

Power-law correlations in both space and time are at least required in order to verify dynamical analogies 
between different catastrophic events. Hence, one may ask how the fluctuations included in a single EQ 
/ ES precursory signal correlate in time. We investigate a possible temporal clustering focusing on a the 
existence of a common potential power-law distribution of burst lifetime (duration) in the populations of: 
(i) electric pulses included in a single ES, (ii) fracto-electromagnetic pulses rooted in the activation of a 
single fault, and (iii) acoustic pulses in laboratory. The result is positive. 
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We claim that the performed analysis captured common principal laws behind the exciting variety of 
complex phenomena under study. 

2. Data Collection 

2.1. Human and rat EEGs 

Electroencephalograms (EEG) are brain signals which provide insight concerning important character- 
istics of the brain activity and yield clues regarding the underlying neural dynamics. 

Human EEG data offered by Andrzejak et al. [30] have been analysed. Two sets, denoted "A" and 
"E", respectively, each one of them containing 100 single-channel EEG segments of 23.6sec duration, were 
employed for this study. Set "A" is comprised of EEGs of healthy volunteers. Set "E" contains seizure 
activities. The segments fulfil the criterion of stationarity [30]. 

EEG signals available from the on-line PhysioNET database (http : //physionet . org/pii6/chbmit/), 
collected at the Children's Hospital Boston, have also been analyzed. This database consists of EEG 
recordings from pediatric subjects with intractable seizures. More details can be found in [31]. In contrast 
to the data found in [30], where set "E" contain exclusively seizure activities, these human EEG correspond 
to signals acquired throughout the whole evolution of the phenomenon (i.e., capturing the electric activity 
before, during and after the seizure). 

Moreover, adult Sprague-Dawley rats were used to study the epileptic seizures in EEG recordings [32, 33]. 
A bicuculline i.p. injection was used to induce the rat epileptic seizures, while EEG signals were recorded. 
The seizure onset time was determined by visual identification of a clear electrographic discharge, and then 
looking backwards in the record for the earliest EEG changes from baseline associated with the seizure. 
The earliest EEG change is selected as the seizure onset time. Note that this data were mainly analyzed 
because they correspond to a controlled epilepsy, where the temporal stages (injection, pre-ietal, ictal and 
post-ietal), have been well justified for their origin. The exact time of the seizure provocation was also given. 

2.2. Pre-seismic EM emissions 

A question effortlessly arises whether there is an "EM-seismograph" available which can be used to 
monitor the evolution of a single fault activation process, in analogy to the EEG which is used to monitor 
the evolution of a single seizure activation process. Such a signal has been reported [34]. Fracture induced 
MHz-kHz electromagnetic (EM) fields allow a real-time monitoring of damage evolution in materials during 
mechanical loading. Crack propagation is the basic mechanism of material failure. EM emissions are 
produced by opening cracks. The radiated EM precursors arc detectable both at a laboratory and geological 
scale [35, and references therein]. Especially, accumulated laboratory and field data suggest that the kHz 
EM activity, which is emerged in the tail of the observed EM radiation, is rooted in the fracture of strong 
entities distributed along the beleaguered by stresses fault sustaining the system [36-39]. 

Herein, we mainly refer to the well documented seismogenic kHz EM activities associated with the: 
Athens (Greece) EQ (A/ = 5.9) that occurred on September 7, 1999 [35-38, 40-42]. 

3. Fundamentals of nonextensive statistical mechanics 

Before presenting our arguments for the problem under study, we briefly introduce the nonextensive 
formalism to which we are referring to. If the physical system involves long-range interactions or long-range 
microscopic memory or (multi-)fractal boundary conditions, it can exhibit a quite anomalous thermody- 
namic behaviour, which might even be untraceable within Boltzmann-Gibbs (BG) statistical mechanics. To 
overcome at least some of these pathological situations, an entropic form Sq has been proposed by Tsallis 
[1, 43], which yields a generalization of standard statistical mechanics and thermodynamics. This entropy 
is defined as: 




(1) 
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where, pi are the probabihties associated with the microscopic configurations, W is their total number, 
and k IS a, positive constant (from now on taken to be unity, without loss of generality). The value 5 = 1 
corresponds to the standard, extensive, BG statistics. 

A property which characterizes the above generalized entropic form is the following: if we have two 
independent systems A and B such that pf^ = pf + pf , the Tsallis entropy satisfies 

Sg{A + B) = Sg{A) + Sg{B) + (1 - q)Sg{A)Sg{B). (2) 

The last term on the right side of the last equation brings the origin of nonextensivity of the resulting 
generalized statistical mechanics. The absolute value of 1 — g indicates the degree of nonextensivity in a 
complex system. If g > 1, = 1 or < 1, 5'^ is subextensive, extensive or superextensive, respectively. The 
index q appears to characterise universality classes of nonadditivity, by phrasing this concept to what is 
done in the standard theory of critical phenomena. Within each class, one expects to find infinitely many 
dynamical systems [1]. 

Perhaps two of the most sound and rich examples of the dynamics of a complex system in crisis are the 
behaviour of brain / earth crust during the ES / EQ generation. A central property of their generation is the 
occurrence of large-scale collective behaviour with a very rich structure, resulting from repeated nonlinear 
interactions among the constituents, namely, firing neurons / opening cracks, of the system. Consequently, 
the nonextensive statistical mechanics [1] is the appropriate framework in order to investigate the process 
of launch of the two shocks under study. 

4. Distinguishing the extreme events in terms of Tsallis entropy 

A basic step in the study of natural extreme events is to clearly distinguish them in the associated time 
series. The Tsallis entropy clearly discriminate the emergence of shocks in EEGs and kHz EM-seismograms 
by means of order of organization. Indeed, a way to examine transient phenomena in a time-series is to 
analyze it into a sequence of distinct time windows of short duration and compute various measures of the 
degree of organisation content in each one of them. We employ here the Tsallis entropy in its symbolic 
form for word length 5 and q = 1.7 (see ref. [42]). A significant Tsallis entropy "drop" accompanies the 
appearance of two extreme events under study. More precisely: 

Fig. 1 shows the EEC of a rat along with its analysis in terms of Tsallis entropy. Three distinct phases 
in this time series were identified. The first (black) phase is the normal (healthy) state, before the applied 
injection which evokes the ES. It is followed by the pre-ictal (green) state which has higher organization 
in comparison to the organization of the healthy state. The first stage of the third (red) phase (ES) has a 
much higher degree of organization even in respect to those reported in the second phase. It is noted that 
the last stage includes the return to the normal state, while the Tsallis entropy gradually increases up to 
that of the healthy state. 

Fig. 2 refers to the kHz EM emission associated with the Athens EQ. Three distinct phases of time series 
were also identified. The first (black) part corresponds to the normal state (EM background) , it is far from 
the EQ occurrence. It is followed by an epoch including a population of EM events sparsely distributed in 
time, which have higher organization in comparison to the organization of the EM background. The launch 
of the third epoch from the background is characterized by an ensemble of intense EM events, which have 
a much higher degree of organization even in respect to those reported in the second epoch. The last EM 
emission (comprised of two bursts) is thought to be due to the fracture of the family of large high-strength 
entities distributed along the fault sustaining the system [35-38, 40-42]. 

Note that both the examined extreme events (red epochs in Figs. 1 and 2), have also been characterized 
by persistency, while the normal states (black epochs) have been characterized by antipersistency [33] . The 
appearance of a high organization dynamics, which is simultaneously characterised by a positive feedback 
mechanism indicating a strong infiuence of excitation of an event on succeeding events, is consistent with the 
emergence of a catastrophic phenomenon and clearly discriminate the pathological from the healthy state. 
This suggestion is further verified by the analysis of human EEGs. Indeed, Fig. 3 shows a human EEG 
along with its analysis in terms of Tsallis entropy and the Rescaled Range analysis (R/S), namely the Hurst 
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Figure 1: The top graph depicts a rat EEG containing an evoked seizure. The arrow shows the time of 
injection. The black, green and red parts refer to the normal, pre-ictal and seizure epochs, correspondingly. 
The last stage includes the return to the normal state. The bottom graph shows the temporal evolution of 
Tsallis entropy, at non-overlapping windows of 1024 samples. 

exponent [44, 45]. The red part corresponds to the ictal phase, as has been annotated in the PhysioNET 
database [31]. It is observed that the aforementioned scheme is also verified here. Characteristically, analysis 
in terms of Tsallis entropy, shows that the seizure part of the signal has a higher degree of organization than 
the rest of the signal. Additionally, as concerns the analysis in terms of the Hurst exponent, the seizure part 
presents persistent behaviour (H > 0.5) in contrast to the rest of the signal. 



5 



2000 



1000 







0.6 



S 0.4 



0.2 





1 1 1 






• • • »* : 
• 


..• • • y: 

* 

«•. • 
1 • 


f » •? 

. *• ' ». 



500 



1000 



1500 2000 
# of windows 



2500 



3000 



Figure 2: The EM time series associated with the Athens EQ recorded by the 10 kHz magnetic sensor. The 
black part refers to the background noise, while the green and red parts refer to the two distinct epochs 
of the emerged seismic EM activity (see text). The bottom graph shows the temporal evolution of Tsallis 
entropy, using non-overlapping windows of 1024 samples. 
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Figure 3: The top graph depicts a human EEG that contains a seizure (red color), as given from the Phys- 
ioNET database [31]. The bottom graph, shows the temporal evolution of Tsallis entropy using successive 
non-overlapping windows of 1024 samples each. 

The higher degree of organization of the patients EEG compared to that of the healthy ones is further 
verified by analysing human EEGs. In Fig. 4 we present two healthy (black) and two patients (red) EEGs 
included in the set A and E, respectively [30]. 

Fig. 5, shows the Tsallis entropies for 100 healthy and 100 patient human EEGs included in the sets 
A and E [30]. As expected, our results depend upon the Tsallis q value. Figure 5, clearly illustrates the 
superiority of the q values restricted in the range 1 < g < 2 to magnify differences of the Sq. It is worth 
mentioning that the nonextensive q parameter values which clearly quantify the temporal evolution of the 
complexity in the EEGs time series are in full agreement with the upper limit q < 2 obtained from several 
studies involving the Tsallis non-extensive framework [46, and references therein]. Moreover, they are in 
harmony with an underlying sub-extensive system, q > 1, verifying the emergence of strong interactions 
in the brain, especially during the occurrence of an ES. However, it is expected that, for every specific 
system, better discrimination will be achieved with appropriates ranges of q values [43]. Thus, a challenge 
will be to estimate the appropriate value of q which is associated with the generation of ESs. We attempt 
an estimation of the appropriate choice of the q index in Fig. 5. We find that the appropriate q index has 
a value of g ~ 1.55. 
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Figure 4: The top graph, depiets a sequence of two healthy human EEGs and a sequence of two human 
epileptic seizures. The healthy EEGs are black coloured (S1,S2) whereas the patient EEGs are red coloured 
(S3,S4). 
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Figure 5: Shannon entropy (top-left), and Tsallis entropy of 100 healthy and 100 patient EEGs belonging 
to the "A" and "B" set, correspondingly [30], for different values of q. The healthy EEGs are black colored 
while the patient EEGs are red coloured. 
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5. Dynamical analogy in terms of nonextensive earthquake dynamics model 

A model for EQ dynamics has been recently introduced by Sotolongo-Costa and Posadas [15]. The 
mechanism of triggering EQs is established through the combination of the irregularities of the fault planes, 
on one hand, and the fragments between them, on the other hand. The fragments size distribution function 
comes from a nonextensive Tsallis formulation, starting from first principles, i.e., a nonextensive formulation 
of the maximum entropy principle. This nonextensive approach leads to a Gutenberg-Richter type law for 
the magnitude distribution of EQs. Silva et al. [29] have subsequently revised this model considering the 
current definition of the mean value, i.e., the so-called g-expectation value as follows, 

log [N (> M)] = log TV + (^—f ) 

where, N is the total number of EQs, N{> M) the number of EQs with magnitude larger than Ad , and 
M K, loge. a is the constant of proportionality between the EQ energy, e, and the size of fragment. 

Now, we examine whether Eq. (3) describes sequences of: different EQs occurred in different faults in 
various seismic regions, fracto-EM events associated with the activation of a single fault, and electric events 
included in a single ES. 

(i) Eq. (3) provides an excellent fit to different regional scismicities each of them including a population 
of earthquakes originated in different faults of the region under consideration. Importantly, the g-values 
are distributed within the range of (1.6 — 1.8) in various seismic regions studied [15, 29, 47, 48]. Recently 
Telesca [17] reported a Maximum Likelihood estimation of the nonextensive parameters of the EQ cumulative 
magnitude distribution. 

(ii) The nonextensive formula (3) describes the detected EM precursors associated with the activation 
of a single fault [38, 40]. We briefly focus on this point. 

The notion of "EM-EQs": Supposing that the background (noise) level of the EM time series A{ti) 
is Anoise, we regard as amplitude Af^m of a candidate "fracto-EM emission" the difference Afemiti) = 
A{ti) — Anoise- We consider that a sequence of k successively emerged "fracto-EM emissions" Afem{ti), 
i = 1, . . . , k represents the EM energy released, e, during the damage of a fragment. We shall refer to this 
as an "electromagnetic earthquake" (EM-EQ). Since the squared amplitude of the fracto-EM emissions is 
proportional to their power, the magnitude M of the candidate EM-EQ is given by the relation: 

M = loge ~ log (J2 [AfemiUf) , (4) 

It has been shown that Eq. (3) fits to the sequence of "EM-EQs" associated with the activation of 
a single fault, for example, Athens and L'Aquila EQs [38, 42, 49, 50]. Herein, iV(> M) the number of 
"EM-EQs" with magnitude larger than M, and a the constant of proportionality between the EM energy 
released and the size of fragment. We note that, the best- fit parameter for these two precursors is given 
by q ^ 1.80, which is clearly in harmony with the exponents obtained in studies of regional scismicities 
[15, 29, 38, 48, 51, 52]. This result means that the activation of a single fault behaves as a "reduced 
image" of the regional seismicity as it was expected. Indeed, From the early work of Mandelbrot [53] , the 
aspect of self-affinc nature of faulting and fracture is widely documented from field observations, laboratory 
experiments, and studies of failure precursors on the small (laboratory) and large (EQ) scale. 

(iii) We concentrate now on single human ESs. Figs. 6a, 6b and 6c show that the distribution of 
magnitudes of electric pulses included in three single ESs available from PhysioNET database ^ [31], follow 
the nonextensive Eq. (3) with q = 1.70 ± 0.03, q = 1.77 ± 0.02 and q = 1.76 ± 0.01, respectively Note that 
G (> M) vs M is depicted in Figs. 6a-6c for convenience, where G (> M) ~ N {> M)/Ntot. , with Ntot. the 
total number of events. 



^ http://physioiiet.org/pn6/chbmit/chb01/chb01_04.6df, Channel 14; F8-T8 / 
http://physionet.org/pn6/chbmit/chb01/chb01_04.edf, Channel 22: FT10-T8 / 
http://physionet.org/pn6/chbmit/chb03/chb03_03.edf, Channel 14: F8-T8. 
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For reasons of completeness and in order to meet the analysis performed by Osorio et al. [28], we also 
refer to the traditional Gutenberg-Richter law (G-R) which states that the number of EQs with magnitude 
greater than M follows the relation: 



log {N> M)=a- bM, 



(5) 



where a and b are constants describing the regional level of seismicity and relative size distribution of events 
respectively. 

We note that Sarlis et al. [16] have shown that above some magnitude threshold Eq. 3 reduces to the 
G-R law with 

--(^) («) 

The experimental data shown in Fig. 6c were also directly fitted to G-R law for events with magnitude 
larger than 4.2, yielding b = 0.63 ± 0.01 (Fig. 6d). Note that the g— value resulting from by fitting Eq. (3) 
to the specific data (Fig. 6c, q = 1.76 ± 0.01) leads to an estimated, by Eq. (6), best = 0.63 ± 0.03 which 
coincides with that deduced via the direct use of the G-R law. 
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Figure 6: (a)-(c) We used Eq. (3) to fit the distribution of the magnitudes of electric events included in three 
single human epileptic seizures available from PhysioNET database [31]. (d) Fitting of the data analyzed in 
Fig. 6c by the G-R law. 
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Figs. 6a-6d show a superiority of the nonextensive Eq. (3) to fit the data compared to the G-R law. 
There is always a elear drop off away from the linear trend of G-R law toward smaller magnitude values 
either we refer to natural or laboratory seismicity; this is also observed for the analyzed ES data. On the 
contrary, the nonextensive Eq. (3), as can be seen from Figs. 6a-6c, is capable of describing the magnitude 
distribution of electric pulses included in a single ES in the whole magnitude range. We emphasize that the 
aforementioned advantage of Eq. (3) in respect to the G-R law is also observed in the description of natural 
seismicity [15, 29]. Hence, Eq. (3) provides the same physics for all the scales. 

Now we refer to the human ESs included in the "E" ^ set [30]. The analysis of the 100 ESs of set 
E reveals that the vast majority of them (84% of the analyzed ESs) can be fitted by Eq. (3) yielding 
q e [1.27, 1.67] (Fig 7f), while the rest of them seem to follow quite different statistics than the one described 
by cither Eq. (3) or the G-R law (e.g.. Fig. 7e). Among the successfully fitted cases, two different groups 
should be discriminated. In the first one, the distribution of electrical pulses included in a single ES is 
nicely described by the Eq. (3) at all scales; Fig. 7a shows such an example. In the second group, the 
data can not be described by Eq. (3) at all scales; the fit is poor lower than a magnitude threshold value 
(refer to the deviation that the experimental data show from the estimated curve in the "knee" region in 
Fig. 7c). However, attempting a series of fittings by Eq. (3) using a different threshold values each time, 
always leads to the same estimated (j— value, although with different accuracies (Fig. 7c). Note that the 
g— values resulting from the fitting of Eq. (3) on the experimental data (for both Fig. 7a and Fig. 7c we 
obtain q = 1.67 ± 0.02) lead to estimated 6— values, by Eq. (6), 6est = 0.99 ± 0.09. The close agreement of 
the estimated and directly calculated 6— values for the two different ES examples analyzed, h = 1.03 ± 0.01 
(Fig. 7b) and h = 0.92 ± 0.01 (Fig. 7d), respectively, verifies that the nonextensive model of Eq. (3) can 
successfully describe the analyzed ES data. 

The observed differentiation of the behaviour of SEs belonging to the revealed two different groups is 
not unexpected. It has been emphasized the great variability in the morphology and mode occurrence of 
seizures, and the great variability of non-ictal EEG patterns [54]. Note that the 100 ESs of set E come 
from different patients and different recording regions of the brain exhibiting ictal activity. Therefore, the 
observed differentiation in fitting Eq. (3) may be a refiection of the variability in the morphology of the 
analyzed ESs. Note that this hypothesis is further enhanced by the observed variety of organization among 
the different signals of set E (Fig. 5). Moreover, although the level of experimental noise is not given 
for the signals of set E, it is expected that different signals (from different patients / electrodes) may be 
characterized by different signal to noise ratios (SNRs). For the signals with lower SNR it is likely that the 
validity of events of magnitude lower than a specific threshold is questionable. 



^http : //eplleptologie-bonn. de/cms/upload/workgroup/lehnertz/S .zip 
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Figure 7: (a) The distribution of the magnitudes of electric events included in a single human seizure included 
in the set "E" [30] showing perfect agreement to the model expressed by Eq. (3), (b) Fitting of the data 
analyzed in Fig. 7a by the G-R law. (c) Experimental data from a single human seizure included in the set 
"E" showing tolerable agreement to the model expressed by Eq. (3), (d) Fitting of the data analyzed in Fig. 
7c by the G-R law. (e) Experimental data not following either Eq. (3) or G-R law (f) The values obtained 
by each one of the 84 out of 100 seizures belonging to the "E" set which could be fitted by Eq. (3). 
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The above presented results support the suggestion that a dynamical analogy in terms of nonextensive 
statistical mechanics exists between ES and EQ generation. Though intriguing to some extent, the above 
mentioned results suggest that a more exhaustive study of the aforementioned biological and geophysical 
shocks in terms of nonextensive statistics is needed to give a deeper interpretation of their generation. 

6. Dynamical analogy by means of "Gutenberg-Richter scaling law" 

In terms of energy the G-R law states that, the probability density function of having an EQ energy E is 
denoted by the power-law P{E) ^ where _B 1.4 — 1.6 [55]. The probability of an ES in a population 
of different events having energy E is proportional to E~^ , where B ~ 1.5 — 1.7 [28]. 

Herein, we examine whether the sequences of electrical pulses included in a single ES follow a power-law 
P{E) ~ E~^ , with a compatible exponent. Fig. 8a, shows that the energies, E, of the electrical pulses 
included in the single rat seizure depicted in Fig. 1 follows the power-law A''(> E) ~ or equivalently, 

the power-law N{E) ^ E^^'^^. Fig. 8b, depicts that the energies, E, of the electrical pulses included in a 
human ESs follows the power-law iV(> E) ^ or equivalently, the power-law N{E) ^ E~^-'''^ . 




Figure 8: (a) The distribution of energies, E, of the electrical pulses included in the rat seizure depicted in 
Fig. 1 follows the power-law iV(> E) ^ E~^-^^ . (b) The distribution of energies, E, of the electrical pulses 
included in a human ES, depicted in Fig. 3, follows the power-law iV(> E) ^ 

As regards the case of the kHz EM precursor associated with the Athens EQ (Fig. 2, upper panel), it 
has been shown that the cumulative number 7V(> A) of pre-seismic EM pulses having amplitudes larger 
than A follows the power- law A^(> A) ^ ^-0.62 |-ggj_ ^j^^g^ probability of an EM-pulse having energy 
E is proportional to E~^-^^ [57]. 

The above mentioned results indicate that sequences of: (i) fracto-EM-pulses included in a single EM- 
prccursor associated with the activation of a single fault, (ii) electric pulses included in a single ES, (ii) 
different EQs occurred in various areas included many faults, and (iv) different ESs, follow the power-law 
P{E) ^ E~^ with rather compatible B-exponents. Notice, in general, differences in constituting elements 
(organic vs inorganic), in scale, and in other properties between the earth and brain may account for 
dissimilarities in the values of exponents [28]. Thus, the reported dynamical analogy in [28] between ESs 
and EQs by means of energy in populations of different ESs and EQs is extended to smaller scale, namely 
in the activation of a single ES / EQ. 
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Importantly, the aforementioned dynamical analogy is extended up to the scale of laboratory scismicity. 
Acoustic / EM emission in rocks has been studied as a model of natural scismicity. Rabinovitch et al. 
[58] have studied the fractal nature of EM radiation induced by rock fracture. Their analysis of the pre- 
fracture EM time series reveals that the cumulative distribution function of the amplitudes follows the 
power A^(> A) A^^-^^, and, consequently, the distribution function of the energies follows the power-law 
P{E) ^ E~^-^^, as it happens in the case of fracture EM radiation associated with the activation of a 
single natural fault. Petri et al. [59] found a power-law scaling behaviour in the acoustic emission energy 
distribution with B = 1.3±0.1. Houle and Sethna [60] found that the crumpling of paper generates acoustic 
pulses with a power-law distribution in energy with S = 1.3 — 1.6. On the other hand, Cowie et al. (1993) 
[61]; Sornette et al. (1994) [62]; Cowie et al. (1995) [63], have developed a model of self-organized EQs 
occurring on self-organized faults. Their theoretical study suggests that the value should he B = 1.3. 



7. Dynamical analogy in terms of waiting times 

Power-law correlations in both space and time are at least required in order to verify dynamical analogies 
between different catastrophic events. Hence, one can ask how the EM / electrical fluctuations included 
in a single EM precursor / ES correlate in time. We investigate the aforementioned temporal clustering in 
terms of burst lifetime (duration) focusing on a potential power-law distribution. We note that in [28], the 
probability-density-function for intervened intervals r was calculated for a population of different ESs. The 
statistic follows a power-law distribution ^ l/(r^+''), with /3 ^ 0.5 for interseizure intervals. 

Figs. 9a and 9b, depict the distributions of the lifetimes of the electric pulses included in a single rat 
(shown in Fig. 1) and human ES (shown in Fig. 3), correspondingly. We observe that both time-series 
follow a power-law distribution ~ I/tiu^ '' / ~ I/t^^-^, respectively. 

We refer to the pre-seismic EM activity associated with the Athens EQ (sec Fig. 2 upper panel). The 
analysis reveals that the durations (lifetimes) display a power-law distribution ~ I/t^,^'^ (Fig. 9c). 
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Figure 9: (a and b) Distributions of the lifetimes of the electric pulses included in a single rat (shown in Fig. 
1) and human ES (shown in Fig. 3), correspondingly, (c) Distributions of the lifetimes of the EM pulses 
included in the precursory kHz EM activity associated to Athens earthquake. The time durations (livetime) 
of the emerged EM events, follows the power-law distribution ~ 1/t^'^. 

Note that Vespignani et al. (1995) [64] measured an exponent = 1.6 via acoustic signals from laboratory 
samples subjected to an external stress. The above mentioned findings further enhance the view that a 
dynamical analogy exists between seismic and epileptic events. 
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8. Discussion &: Conclusions 



The field of study of complex systems considers that the dynamics of complex systems are founded 
on universal principles that may be used to describe disparate problems ranging from particle physics to 
economies of societies. This is a basic reason for our interest in complexity [2-6, 6-13]. 

Authors have suggested that dynamics of EQs and neurodynamics can be analyzed within similar math- 
ematical frameworks [25-27] . The present paper reports the first indication that a dynamical analogy exists 
between the generation of a single epileptic seizure and a single earthquake activation. The results of this 
work were built on the concepts of nonextensive statistical mechanics. More precisely, a recent nonextensive 
model for earthquake dynamics leads to a Gutenberg- Richter type law for the magnitude distribution of EQs 
(formula 3). We showed that the populations of fracto-EM-pulses rooted in the fracture of the backbone 
of strong entities distributed along a single fault sustaining the system, and electric pulses included in a 
single epileptic seizure follow the above mentioned nonextensive distribution with compatible nonextensive 
g-parameter. 

Our analysis also revealed common "pathological symptoms" of a transition to the generation of seem- 
ingly different two extreme events under study. The generation of a single epileptic seizure/ earthquake is 
accompanied by the appearance of: (i) significant lower non-extensive Tsallis entropy, namely, higher orga- 
nization, and (ii) persistency, i.e., a mechanism which is characterized by a positive feedback, thus leading 
the complex systems out of equilibrium. 

The nonextensive formula (formula 3) also describes the magnitude distribution of solar flares and mag- 
netic storms with well-matched nonextensive g-parametcr. The two last extreme events also are accompanied 
by the appearance of the above mentioned "pathological symptoms" , namely, higher organization and per- 
sistency [33-37, 65-67] . We also investigated the existence of a dynamical analogy between a single epileptic 
seizure / earthquake generation by means of scale- free statistics, namely, the traditional Gutenberg- Richter 
distribution of event sizes and the distribution of the waiting time until the next event. Our study offers 
and justifies a positive answer. One universal footprint seen in complex systems is self-affinity. The results 
reveal that the populations of (a) earthquakes occurred in different faults, (b) epileptic seizures occurred in 
different patients, (c) fracto-EM-pulses rooted in the fracture of the backbone of strong entities distributed 
along a single fault sustaining the system, and (d) electric pulses included in a single epileptic seizure follow 
the same scale-free statistics. Importantly, the analogy is extended up to the laboratory scale of fracture. 

It has reported that electromagnetic flashes of low-energy 7 — rays emitted during multi-fracturing on 
a neutron star, electromagnetic pulses emitted in the laboratory by a disordered material subjected to an 
increasing external load, and electromagnetic pulses emitted in the fleld during the activation of a fault share 
distinctive statistical properties with earthquakes, such as power-law energy distributions [5, 6, 24, 36, 58, 68]. 
The neutron starquakes may release strain energies up to 1046 erg, while, the fractures in laboratory samples 
release strain energies approximately a fraction of an erg. An earthquake fault region can build up strain 
energy up to approximately 1026 erg for the strongest earthquakes. Consequently, the limits of similarity in 
the dynamics of different multi-fracturing extreme events are dramatically expanded. 

In summary, our experimental results support previously reported theoretical and experimental argu- 
ments [25-27] that dynamics of earthquakes and neurodynamics can be analyzed within similar mathematical 
frameworks, while the dynamical analogy by means of non-extensive and traditional free scale laws between 
earthquakes and epileptic seizures is extended up to solar flares, magnetic storms and neutron starquakes 
dynamics. 
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